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ABSTRACT 

Microstructural images may be "tone pulse encoded" and subsequently 
Fourier transformed to determine the two-dimensional density of frequency 
components. A theory is developed relating the density of frequency components 
to the density of length components. The density of length components 
corresponds directly to the actual grain size distribution function from which 
the mean grain shape, size, and orientation can be obtained. 

INTRODUCTION 

Material characteristics, such as tensile strength, hardness, yield 
stress, fracture stress, impact resistance, and ultrasonic attenuation^ 
are directly related to the grain size distribution in polycrystalline 
materials. Thus prediction of these properties requires detailed knowledge of 
the grain size distribution. While the theoretical determination of grain 

D 

size distribution has received considerable attention. The experimentally 
measured distribution function has received limited acceptance, primarily due 
to complexity of real microstructures where the grain topology varies 
considerably. Hence researchers have acquiesced with determination of the 
mean or average grain diameter which Ideally should be determined from the 
grain size distribution. 

O 

There are currently several accepted techniques for determining the 
mean grain size without measuring the grain size distribution function. Host 
methods rely on a line Intercepts along random directions or in circular paths; 



alternatives estimate the circumference or area of each grain. .Current 
techniques are not applicable to any arbitrary system and the researcher must 
often determine which method will yield the most accurate data. The guidelines 
for making this decision are qualitative generalizations and, as such, can lead 

g 

to considerable error in the results. 

The most common method used for determining mean grain size are the ASTM 

g 

line intercept and the ASTM comparison method. The line intercept method 
has been shown to work well for materials having regular grain microstructure; 
however, this method is not applicable to materials exhibiting oriented or 
complex, irregular structure, and when applied by different researchers, the 

g 

line intercept method can lead to results that vary up to 50 percent. The 
comparison method relies on the researchers judgment in comparing standard ASTM 
grain charts with the observed microstructure. Neither of these methods is 
used for determination of the grain size distribution. 

This paper describes a technique for determining the grain size 
distribution function from a metallographically prepared sample. By applying 
two-dimensional Fourier transform theory to tone pulse encoded microstructural 
images, the grain size distribution function is determined. The resulting 
relationship is two-dimensional and yields the grain size distribution, mean 
grain shape, size and orientation. 


THEORY 

One Dimensional Fourier Transforms 
The Fourier transform of an arbitrary function f(x) is given by.^ 

•f-00 

(1) F(u) = f f(x) dx 

—00 

|F(u)| represents a density function such that the number of components 
contained in the interval |du| is given by^^ |F(u)||du|. The variable x may 
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refer to any quantity such as frequency, time, length, etc. For example, 1f x 
Is time then u Is frequency and lF(u)| represents the density of frequency 
components. Given a frequency u the period x' Is’ given by x' = 1/u. 
Throughout this work we will refer to x' as the reciprocal variable of u 
and u as the reciprocal variable of x'. 

Equation (1) may be expressed with respect to the period variable, x'. 

(2) |P(x')|ldx'| = |F(u)||du| 

where |P(x')|ldx'| Is the number of period components contained In a 
period Interval |dx'l- Then 

(3) |P(x')l = |F(u)l 
using 


(4) 


dx' 

du 



we obtain, 

(5) |P(x')| = |F(u)|u^ 


TWO-DIMENSIONAL FOURIER TRANSFORM THEORY ' 

A two-dimensional Fourier transform contains two variables, such as x 

g 

and y and Is given by 

(6) F(u,v) = [ [ f(x,y) dx dy 

-00 -00 


Here the number of components per area Interval |du||dv| Is given by 

lF(u,v) I |dul |dvl . Equation (6) may also be expressed In polar coordinates. 

1 2 

With dx dy = rde dr we have 


F(s,4>) = f f f(r,e) rde dr 

-'n 


-Lff ^0 

V 2 2 #2 2 

u + V and r = + y 
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The number of components per area Interval 1s 

(8) |F(s,<t.)||sd 4 .||ds| 

and the number of components per radial length at an a-ngle <t> Is 

(9) |F(s, 4 >)l|ds|. 

Expressing Equation 9 In terms of period reciprocal variable r' we have 

(10) |P(r'. 4 >)l|dr'| = |F(s,<t.) | |ds 1 

where |P(r',4>)| Is the density of period components In the Interval |dr'|- 
Using s = 1/r' we obtain 

(11) |P(r'.4.)| = |F(s. 4 >)|s^ 

Determination of Boundary Length Distribution Function 

Suppose there exists a set of barriers placed along the axis as shown In 

Figure 1(a). The boundaries are randomly spaced along the x-axIs and have 

the same amplitudes. We are Interested In determining the distribution of 

lengths between the boundaries over the entire length of the x-axIs. The 

Fourier spectra of a three cycle tone pulse having an amplitude 2A, with x' 

and fundamental frequency u^ = 3/x' Is given by 

! s1n [■"■(u^ - u)x'] s1n[ir(u^ + u)x']1 
[it(Uq - U)X'] “ [it(Uq + u)x'] j 

If a three cycle, sine wave, tone pulse Is formed between adjacent barriers 
(see Figure 1(b)) the resulting waveform will be made up of a series of tone 

I I 

pulses each having fundamental frequencies at 3/x^, where x^ refers to 
the length between the 1 and (1+1) barriers. The Fourier transform of the 
waveform will yield |F(u)|, the density of frequency components, and 
|P(x')| will provide the density of length components (or the length 
distribution function). Three cycles will correspond to a fundamental period 
one third the true width between the barriers (This true period or length will 


(12) |F(u)| = Ax' 
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be referred to as the barrier width.). This 1s easily corrected by rescaling 
the fundamental periods by a factor of three. 

It Is possible for the relative phases between the tone pulses to Interact 
by constructive and destructive Interference. This would lead to an erroneous 
density of frequency components where constructive or destructive phase 
Interference dominates. If the phase Interference Is Incoherent (I.e., random) 
or negligible, then the density of frequency components Is accurate. 

Here It Is assumed that phase Interaction between tone pulses Is either 
negligible or Incoherent. This Is a relatively strong assumption, but It will 
be shown later to be a reasonable one. It Is pointed out that for this work 
we will use a digital 512 point fast Fourier transform (FFT) which limits the 
resolution of the x-axIs to 512 points. 

In order to form a smooth barrier width distribution function 512 random 
barrier systems similar to that shown In Figure 1(a) will be evaluated. .Each 
of the 512 barrier systems Is frequency coded with tone pulses and subsequently 
Fourier transformed to determine |F^(u)| where the subscript n refers 
to the n^^ barrier system. The total density of frequency components for 
the entire set of barrier systems Is given by 

512 

(13) IF (u)| = Z IF (u)| 

^ n=l " 

from which the total density of length components Is obtained 

(14) |Pt(x')| = |F^(u)|u^ 

Figures 2(a) and (c) show the density of length components and the actual 

barrier length distribution function, respectively. For lengths less than 

2 

x' < 13 pixels (u > 0.076), |F(u)|u deviates considerably from the actual 

barrier distribution function. This deviation originates from the 

nonnegllgible higher order harmonics that are present for a three cycle tone 

2 

pulse. These high frequency components are enhanced by the factor u 
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(Eq. (14)). Equation (12) Indicates that the density of frequency components 
at the fundamental frequency of a three cycle tone pulse Is proportional to 
the tone pulse width, x'. If the amplitude of each three cycle tone pulse Is 
diminished In proportion to Its corresponding width, then the resultant density 
of frequency components, at the fundamental frequency for each tone pulse. Is 
normalized to a constant value, A. While the higher order harmonics are 
correspondingly diminished In amplitude. For example, the waveform shown In 

Figure 1(b) Is modified by decreasing the amplitude of each tone pulse by 

' ' th 

1/x^ where x^ Is the width of the 1 tone pulse. If the phases 

of the tone pulses exhibit negligible or Incoherent Interaction then the 

Fourier transform of the modified waveform for the barrier system Is 


(15) iFj(u)i 

1 

where ^ denotes the Fourier transform and f^(x) describes the wavefofm 
of the 1^^ tone pulse. Since x^ Is constant for the 1^*^ tone 
pulse we have 



(16) 





for the barrier system. 

The density of frequency components for a set of 512 barrier systems becomes 


(17) 


|F„(t.) 


J J 1 




l/(f^j(x))| « |F^(U)|U 


where the subscript m refers to the modified system. Using Equation (14) the 
density of length components Is 


(18) 


|P(x')| = |F^(u)|u = |F^(u)|u^ 
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Equation (18) Indicates that when the Fourier spectra of the modified system 
lFjj^(u)| Is multiplied by u the density of length components Is obtained 
(compare Equations (18) and (14)). The results follow similarly In two 
dimensions (from Equation (11)) to yield the density of radial length 
components at an angle 4 ., 

(19) |P(r',<t>)| = |FJs,4.)1s = 1F^(s,4.)|s^ 

The approximation In Equation (17) become equivalence when F(f^(x)) 

approach delta functlons^^ at the fundamental frequencies, u^ = 3/x^. 

In other words, If the number of cycles In each tone pulse Is Increased the 

|F (u)lu, the density of length components, becomes a more accurate 
m 

representation of the true density of length components. Unfortunately, 
Increasing the number of cycles of each tone pulse In a digital 512 point array 
prohibits tone encoding of small lengths. Here, the Nyquist frequency^*^ Is 
the limiting factor. Therefore, there Is a trade off between the required 
measurement of small lengths and the required accuracy of the length 
distribution function. 

The resultant modified waveform Is shown In Figure 1(c). This 

modification Is done for the entire 512 barrier systems. The density of length 

components for the modified barrier set Is shown In Figure 2(b) and Is given 

by |P(x')| . |FJ.)|u where one addUlonal factor of o 1s pre-1nclyded 1n 

the construction of the modified barrier systems (see Equation (18)). This 

2 

modification reduces the high frequency enhancement from u to u and 
yields a more accurate representation of the length distribution function as 
shown In Figure 2(b). The flat region of the distribution for pixel periods 
less than 10 pixels Is Identified as residual noise resulting from higher 
frequency harmonics. This flat noise region Is observed In all systems 
Investigated and Is Intrinsic In the Fourier transform theory. This easily 


7 



Identifiable region Is clipped to a zero value to eliminate the higher 
harmonic noise. Figure 2(d) Is an overlay of curves shown In Figure 2(a), 

(b), and (c). 

The results shown In Figure 2 Indicate that the density of length 
components corresponds directly to the actual length distribution function. 
This supports the assumptions that: (1) The phase Interaction of tone pulses 

Is Incoherent or negligible. (2) The selection of three cycles per length Is 
sufficient to yield an accurate result. 

Application of Fourier Transform Theory to Hicrostructural Images 

To apply Fourier transform theory to microstructural Images we tone pulse 
encode the two-dimensional Image. The construction of a two-dimensional tone 
pulse encoded representation of an arbitrary microstructural Image requires 
enhancement of the original Image as Indicated below. The microstructure, of 
stock nickel 200 Is shown In Figure 3(a). The Image was digitally recorded 
Into a 512 by 512 pixel array via a vidicon camera connected to a video 
digitizer. The Image Is a numerical average of .five digitizations. The gray 
scale resolution of the video digitizer allows for pixel Intensities to vary 
from 0 to 255, where an Image Intensity of 0 and 255 correspond to black and 
white, respectively. The grain boundaries are clearly visible as 
Interconnecting dark curves. 

The grain boundaries are enhanced by taking the two-dimensional digital 

Q 

gradient (Figure 3(b)) of the microstructure shown in Figure 3(a). A gray 
level line scan along the x-directlon (In Figure 3(b)) at y = 255 Is shown 
In Figure 3(c). A series of peaks are observed In this line scan that 
correspond to the position of the grain boundaries. The peaks vary in 
amplitude where the magnitude derivative between grains varies In agreement 
with the optical Intensity. Between the peaks there exist noise generated by 
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variations of the optical Intensity within the grains. This noise Is also 
observed In the original microstructural Image (Figure 3(a)). The within grain 
noise Is removed by determining the magnitude of the derivative of the weakest 
(that Is, the most difficult to determine optically) grain boundary. All 
pixels In the gradient image having a value lower than this value are set to 
zero while all others are set to a constant nonzero value. As a result we 
have a clean Image of the grain structure where each grain boundary Is given 
equal weight (I.e., the same Intensity). A gray level line scan along the 
x-directlon at y = 255 for the noise reduced, grain boundary enhanced Image 
(Figure 4(a)) Is shown In Figure 3(d). The line scan reveals a series of 
equal amplitude, half cycle pulses. Each pulse corresponds to Its respective 
grain boundary along the scan direction. The enhanced Image Is a 
two-dimensional representation of the one-dimensional case shown In Figure 
1(a). In order to assure that there Is equal weighting of grains In all 
directions the enhanced Image Is circularly masked as shown In Figure 4(a). 

Next a tone pulse encoded image containing the fundamental harmonics of each 
grain 1s generated. This Is done as In the one-dimensional case. The tone 
pulse begins at a grain boundary and ends at a grain boundary and has a width 
equal to the width of the grain. A new tone pulse Is started at an adjacent 
grain boundary which generally has a different width. This Is done digitally 
starting at the center of the Image and proceeds outward along a fixed number 
of radii to the perimeter of the circular mask. The process Is repeated for 
all angles (0 to 2ir rad) so that the frequency coded Image appears as shown In 
Figure 4(b). Figure 4(c) 1s a gray level line scan along the x-d1rect1on of 
Figure 4(b) at y = 255 and Is similar to that shown for the one-dimensional 
case ( see Figure 1(c)). 

Nine different photomicrographs of the same sample are tone pulse encoded 
and then digitally Fourier transformed using a complex 512 point, hardwired, 
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fast Fourier transform array processor. The density of length components at 
an angle 4> 1s given by Equation (19) and Is shown In Figure 5(a). This 
figure represents the density of length component as i function of the 
reciprocal length, s. We may better understand this result by plotting the 
result In Figure 5(a) as a function of 1/s (I.e., length) as shown In 
Figure 5(b). Figure 5(b) represents the density of length components as a 
function of length, r', and corresponds directly to the actual grain size 
distribution function |P(r',4>)l- From this figure the grain size distribution 
and mean grain size of the original microstructure may be determined along any 
direction. 

The mean grain length 0(4>) along an arbitrary angle 4> Is obtained from 
Figure 5(b) and Is given as 


( 20 ) 


ElP(r',4>)| r' 
- ElP(r'.<J>)l 

r ' 


where r' = V? +y^ 

The mean grain shape, size, and orientation obtained from Equation (20) 
is shown In Figure 5(b) (Inset). For this microstructure the mean grain shape 
Is roughly circular with no net orientation. 

The grain size distribution along the x-d1rect1on (4> = 0), |P(r',0)|, 

Is shown In Figure 5(c) as a dark curve. Also shown for comparison Is a 
histogram of the grain size distribution along the x-directlon obtained by 
using a commercially available Image analyzer. The histogram exhibits a slight 
shift In the distribution toward larger grains. This Is most probably due to 
the technique used for determining grain size for this Instrument. The 
histogram obtained Is determined by measuring the maximum length (In the x 
direction) of each grain. Therefore It Is expected that the histogram be 
shifted to toward larger grains. The “tone pulse encoded" and commercial 
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Instrument results yield mean grain diameters (along the x-dlrejctlon) of 60 
and 72 ym, respectively, and differ by 16 percent. 

As the mean grain size of a real microstructure ts not known a priori and 

o 

Is found to be dependent on the methodology used for Its determination , It 
seems appropriate to provide a test system for which the mean grain size, shape 
and orientation are known. A simulated microstructure containing a 
distribution of sizes of Identically shaped and oriented grains shown In 
Figure 6(a) was used. After tone pulse encoding, the space between the grains 
was masked to a zero value tone code to eliminate treating them as "artificial" 
grains. The resultant density of length components Is shown In Figure 6(b). 

The mean grain size, shape and orientation are obtained by use of 
Equation (20). The reconstructed mean grain shape and orientation Is shown In 
Figure 6(c) along with the actual mean grain. The reconstructed grain Is 
properly oriented and of similar shape to the actual mean grain. The largest 
error In the reconstruction of the mean grain has occurred at the regions of 
sharp radii of curvature (top and bottom of the actual grains). The error at 
these regions Is 12 percent. 

This work may also be applied to ceramics for the determination of pore 
size distribution. Figure 7(a) Is a digitally enhanced photomicrograph of a 
sintered ceramic showing a distribution of pores. The area between the pores 
was set to a zero value tone code after tone pulse encoding. The pore size 
distribution function shown In Figure 7(b) exhibits an anisotropic 
distribution. The mean pore (see Figure 7(b) Inset) Is oval In shape and Is 
oriented along the vertical axis with an aspect ratio of 10:7. 

DISCUSSION 

The entire process of digitizing Images, Image enhancement, tone pulse 
encoding, Fourier transformation and determination of the grain size 
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distribution function takes about 15 hr using a nonvirtual minicomputer, 
hardwired array processor (for FFT) and a video digitizer. The operator needs 
to Interact during the digitization of the microstrucfural Images which takes 
about 20 min. Host of the computer time Is spent reading and writing to the 
video digitizer. This time could be easily reduced by at least an order of 
magnitude by using a virtual minicomputer that can perform mathematical 
operations on an entire Image contained In a random access memory array. By 
comparison the commercial Image analyzer requires about 2 hr of operator 
Intensive activity (tracing grain boundaries) to yield grain distribution data 
along two perpendicular directions. 

Since each microstructural Image contains a different set of grains the 
grain size distribution will not be Identical for each Image. Thus averaging 
over a number of photomicrographs^^ Is necessary to estimate the real grain 
size distribution. The Increase In accuracy of the grain size distribution 
function Is proportional to where N Is the number of averaged Images. 

Nine Images provide for an Increase In accuracy by a factor of three; It 
should be noted that It would take 100 Images to Increase the accuracy by a 
factor of 10. 

In order to tone pulse encode three full cycles between adjacent 
boundaries, at lease 10 pixels are required. Therefore, the smallest grain 
width measurable has a length of 10 pixels, and all grains having widths less 
than 10 pixels are not Included In the determination of the grain size 
distribution. This Is observed In Figure 5(c) where the grain size 
distribution begins a sharp rise at 10 pixels. The functional form of the 
grain size distribution Is not known for smaller grain widths. It could be 
argued that more data would be available by tone pulse encoding with one or two 
cycles Instead of three cycles. By using one or two cycles In the tone pulse, 
grains having widths of 4 or 7 pixels, respectively, could be encoded. In 
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order to address this argument we must examine the gradient Imag.e (Figure 4(a)) 
where the uncertainty of the grain widths Is 3 pixels. 

Tone pulse encoding with fewer cycles would allow for the Inclusion of 
smaller width grains (grains having widths 4 to 7 pixels). However, the 
current uncertainty (3 pixels) In the grain width would remain the same making 
detailed measurement of these narrow width grains of questionable validity. 

The percent uncertainty In the grain width Is given by 

X 

l^xl = ^ • 100 percent, where Ax = 3 pixels and x Is the grain width. 

The percent uncertainty Is Inversely proportional to the grain width. That Is, 
the larger the grain width to be measured the more accurate the measurement. 
There Is less that 15 percent uncertainty (due to the uncertainty In the grain 
widths) In the grain size distribution function for grain widths greater than 
20 pixels. 

This work could be enhanced In several ways. A hardwired tone pulse 
encoding processor could be used In conjunction with a virtual 
superminicomputer to reduce the processing speed to a few minutes. Increased 
accuracy (via averaging) of the distribution function could be available with 
a computer controlled vidicon camera scanning the actual microstructural 
surface. Also, an Image array greater than 512 by 512 pixels would allow for 
finer definition of boundaries allowing for tone encoding of narrower widths. 

This technique Is not limited to grain or porosity distributions. It may 
also be applied to any system for which the length distribution function Is 
required. Some specific examples are the determination of whisker size 
distribution function In chopped fiber composites as well as the length 
distribution function between whiskers. This also applies to continuous fiber 
composites. In metals the Inclusion and metallic phase length distribution 
functions may be determined. Application of this work to three mutually 
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perpendicular planar sections will yield the full three-dimensional mean grain 
and/or pore size, shape, and orientation. 

CONCLUSIONS 

Microstructural Images may be tone pulse encoded and subsequently Fourier 
transformed In order to determine the two-dimensional density of frequency 
components. A theory has been developed and applied to relate the density of 
frequency components to the density of length components. The density of 
length components corresponds directly to the two-dimensional grain and/or pore 
size distribution function from which the mean grain or pore size, shape, and 
orientation are obtained. 
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